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^ ' Abstract 

O' 

•^ ' In this article, we assume that there exist hidden charmed tetraquark 

^^ . states with the spin-parity J^ = 1~ , and calculate their masses with the QCD 

I sum rules. The numerical result indicates that the masses of the vector hidden 

<~~^- charmed tetraquark states are about Mz = (5.12ib0.15) GeV or Mz = (5.16± 

o . 0.16) GeV, which are inconsistent with the experimental data on the tt'^XcI 

' I invariant mass distribution. The hidden charmed mesons Zi, Z^ or Z may 

Q ' be scalar hidden charmed tetraquark states, hadro-charmonium resonances or 

|-G^. molecular states. 
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^. ■ 1 Introduction 

c^ \ The Babar, Belle, CLEO, DO, CDF and FOCUS collaborations have discovered (or 

O ■ confirmed) a large number of charmonium-like states [IlElElllllHlEllTllHllHl [lOl [TT] , 

K»- ■ and revitalized the interest in the spectroscopy of the charmonium states p^ [T3l 

><: naE]. 

H \ The X(3940) decaying into D*D and the Z(3930) decaying into DD have been 

■ - - ' tentatively identified as candidates for the missing charmonium states t]"^ and Xc2 

respectively. The X(3872) decaying into 'k+tx~J/^Ij, 7r+7r-7r°J/^, the F(4260), 
y(4008) decaying into tt+tt-J/^, the r(3940) decaying into uJ/^, the F(4325), 
y(4360), y(4660) decaying into 'k^'k~%1)' have odd properties comparing with the 
expectations of the charmonium models [121 [13 O [15] . 

Many possible assignments for those states have been suggested, such as multi- 
quark states (whether the molecular type [HI [T7] or the diquark-antidiquark type 
[iHl [ini [20] ) , hybrid states [211 [221 [23], charmonium states modified by nearby thresh- 
olds [211 [25], threshold cusps [26], etc. 

The observed decay channels are J/tp7r^n~ or ijj'7r~^7r~, an essential ingredient for 
understanding the structures of those mesons is whether or not the vrvr comes from 
a resonance state. For example, there is an indication that the y(4660) has a well 
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defined intermediate state in the vrvr invariant mass distribution, which is consistent 
with the scalar meson /o(980) [27], the F(4660) can be taken as a fo{980)ip' bound 
state [28j; though other interpretations such as a baryonium state [29j or a canonical 
5^Si cc state [30] are not excluded. 

The ^+(4430) observed in the V^V"*" decay mode is the most interesting subject 
[31]. It can't be a pure cc state due to the positive charge. There are many theoretical 
interpretations for its structures, such as the hadro-charmonium resonance [HI [32] , 
the molecular D*Di{D[) state [531 [Ml [2S1 [Ml [SZl [3H] , the tetraquark state [311 iDl 
HH HSl Us], the cusp in the D*Di channel [H], the radially excited state of the Dg 
[l5], etc. We can distinguish the multiquark states from the hybrids or charmonia 
with the criterion of non-zero charge. 

Recently the Belle collaboration reported the first observation of two resonance- 
like structures (thereafter we will denote them as Zi and Z2 respectively) in the 
ir'^Xci invariant mass distribution near 4.1 GeV in the exclusive B^ -^ K^tt^Xci 
decays [^. Their quark contents must be some special combinations of the ccud, 
just like the Z~^{A430), they cannot be the conventional mesons. The Breit-Wigner 
masses and widths are about Mi = 4051 ± 141^0 MeV, Ti = 821^1+47 MeV, M2 = 
4248+^^+^f MeV and T2 = 177+^^+^f MeV. The significance of each of the n+Xci 
structures exceeds 5a, including the effects of systematics from various fit models. 

The Z (denote the Zi and Z2) lie about (0.5—0.6) GeV above the vr+Xci threshold, 
the decay Z -^ tt+Xci can take place with the "fall-apart" mechanism and it is OZI 
super-allowed, which can take into account the large width naturally. The spins of 
the Zi and Z2 are not determined yet, they can be scalar or vector states |46j . 

If they are scalar mesons, the decays Z -^ vr+Xci occur through the relative P- 
wave with the phenomenological lagrangian C = gx'^i'^daZ — ZdaT^). On the other 
hand, if they are vector mesons, the decays occur through the relative S'-wave with 
the phenomenological lagrangian C = gx'^Zair. 

The typical decay mode Z -^ D^D^ is kinematically allowed, and the width may 
be comparable with the corresponding ones of the decay mode Z -^ tc^Xci, we can 
determine the spins of the Z with the angular distributions of the final states D^D^ . 
If the decays Z — >■ D^D^ are not observed (or the widths are rather narrow), the Z 
may be hadro-charmonium resonances [11], i.e. bound states of a relatively compact 
charmonium (Xci) inside a light hadron (vr"*") having a larger spatial size, the decays 
Z -^ vr+Xci occur with the "fall-apart" mechanism; the decays Z -^ D^D^ take 
place through the final-state re-scattering effects, Z -^ tt^Xci -^ D^D^ , the widths 
may be very narrow. 

The masses of the D^ and D^ are about Md = 1.87 GeV, the Z may also be P- 
wave D^D^ molecular states with the spin-parity 1^, as the additional contribution 
from the relative P-wave is about 0.5 GeV in the potential quark models, the decays 
Z -^ D^D^ occur through "fall-apart" mechanism and have much larger widths than 
the corresponding decays Z — > D^D^ —>■ 7r~^Xci, which occur through final-state re- 
scattering effects. One may also think that they are D*Di (or D*D[) molecular 
states, Mb. = 2.01 GeV, M^, = 2.42 GeV, M^^ = 2.43 GeV, the bound energy 



is about (0.2 — 0.4) GeV, which may be beyond capacity of the one-vr and one-cr 
exchange. 

The mass is a fundamental parameter in describing a hadron, in order to identify 
the Z\ and Z2 as tetraquark states, we must prove that the masses of the correspond- 
ing tetraquark states he in the region (4.1 — 4.3) GeV. Furthermore, whether or not 
there exist such hidden tetraquark configurations is of great importance itself, be- 
cause it provides a new opportunity for a deeper understanding of low energy QCD. 

In this article, we assume that the hidden charmed mesons Z^ and Z^ are vector 
tetraquark states, which consist of a pseudoscalar (scalar) diquark and an axial- 
vector (vector) antidiquark, and study their masses with the QCD sum rules [ITIHH]. 
As their spins are not determined yet, they may also be scalar hidden charmed 
tetraquark states, we study this possibility with the QCD sum rules j4Si|. 

In the QCD sum rules, operator product expansion is used to expand the time- 
ordered currents into a series of quark and gluon condensates which parameterize the 
long distance properties of the QCD vacuum. Based on current-hadron duality, we 
can obtain copious information about the hadronic parameters at the phenomeno- 
logical side. 

It is difficult to distinguish the mesons Z^ and Z^ with the QCD sum rules 
approach, as the mass gap between them is very small. The mesons Z\ and Z2 lie in 
the region (4.0 — 4.3) GeV, we study whether or not there exist 1~ hidden charmed 
tetraquark states Z in this energy region and possible tetraquark identification of the 
mesons Z^ and Z2. Furthermore, the vr+Xci invariant mass distribution amplitude 
can also be represented by a single Breit-Wigner mass formula, Mz = 4150^'^g MeV 
and r^ = 352+^^ MeV jlS]. 

The article is arranged as follows: we derive the QCD sum rules for the masses 
and pole residues of the Z in section 2; in section 3, numerical results and discussions; 
section 4 is reserved for conclusion. 



2 QCD sum rules for the tetraquark states Z 

In the following, we write down the two-point correlation functions n^j,(p) (denote 
^fiuip) ^^^ ^^uip)) i^ the QCD sum rules, 

Ki'iP) = ^ I d'xe'^-{0\T {j/r^,{x)J/r^l{0)}\0) , (1) 

J^(x) = e*^'=e*'""Mj(x)Ccfc(x)c„(x)7MCrfT(x), (2) 

V^{x) = e^^'=e^™"nJ(x)C75Cfc(a;)c^(x)757MC(iT(a;), (3) 

fzM^e, = {0\J/v,iQ)\Z{p)), (4) 

we choose the vector currents J^{x) {C — 6*7^ type) and rj^{x) (C75 — 6*7^75 type) 
to interpolate the tetraquark states Z, the fz is the pole residue and the e^ is the 
polarization vector. If there exist vector tetraquark states Z in the tt^Xci invariant 



mass distribution, it is convenient to construct the tetraquark currents with the 
pseudoscalar (scalar) diquark and axial-vector (vector) antidiquark, as the mesons 
TT and Xci have the spin-parity 0~ and 1"^ respectively, the decays Z -^ vr+Xci can 
occur through relative S'-wave. We can also interpolate the vector tetraquark states 
with the currents Jj^^x) and rij^{x), 

Jlix) = e'^h'^^uJix)C^M^)cm{x)Cd^{x), (5) 

Vlix) = e*^'^e*™"Mj(x)C7^75Cfc(x)c„(x)75CrfT(^), (6) 

which consist of a pseudoscalar (scalar) antidiquark and an axial-vector (vector) 
diquark. Our analytical results indicate that the current J^{x) (?7^(x)) and J^(x) 
(?7^(x)) lead to the same expression. The special superpositions tJ^{x) + (1 — t) J^^(x) 
and trif^{x) + [1 — t)r]l^{x) cannot improve the predictions remarkably, where t = — 1. 
The correlation functions Il^^{p) can be decomposed as follows, 

ii,Ap) = {-9,. + ^)Mp') + ^'Mp') + ---, (7) 

due to Lorentz covariance. The invariant functions Hi and Ho stand for the contri- 
butions from the vector and scalar mesons, respectively. In this article, we choose 
the tensor structure g^^ — ^^^ to study the masses of the vector mesons. 

Basing on the quark-hadron duality [171 HH] , we can insert a complete series of 
intermediate states with the same quantum numbers as the current operators Jfi{x) 
and ?7^(x) into the correlation functions li^viv) to obtain the hadronic representation. 
After isolating the ground state contribution from the pole terms of the Z, we get 
the following result. 






Ml-p^ 






+ ■■■. (8) 



In the following, we briefly outline operator product expansion for the correlation 
functions n^t,(p) in perturbative QCD theory. The calculations are performed at 
large space- like momentum region p^ -C 0. We write down the "full" propagators 
Sij{x) and Cij{x) of a massive quark in the presence of the vacuum condensates 
firstly HH], 

6' X iS' 'X 

-j^iqgsCrGq) + -^mg{qgscrGq) :/: - 

^G-.(^a- + a-;^) + ..., (9) 



Cij{x) = -nr-rj d ke 



-ik-x 




+ rUc) + {}t + rricjaafi 



{2ttY J H~mc 4 {k'^-m, 

t"^ asGG k'^ + rric^ 



I? 



Tf^ asGG k^ + TTic^ , , 



where {sggCrGs) = {sgscrapG^^s) and ( °°'^'^ ) = ( "" ""^ ), then contract the quark 

fields in the correlation functions Il^vip) with Wick theorem, and obtain the result: 



[X 



Tr [j^CS^,^{-x)C^,C^,U-x)] , (11) 

n^^(p) = i^ijk^imn^i'j'k'^i'm'n' j d'^xc'^'^Tr ^I^C Sl^^,{x)C-i^C jj,{x)\ 

Tr [757^C5'J,„(-a;)C7^75C„'„(-a;)] . (12) 

Substitute the full u, d and c quark propagators into the correlation functions n^j,(p) 
and complete the integral in the coordinate space, then integrate over the variables 
in the momentum space, we can obtain the correlation functions ni(p^) at the level 
of quark-gluon degrees of freedom. Once analytical results are obtained, then we 
can take current-hadron duality below the threshold Sq and perform Borel transform 
with respect to the variable P^ = —p"^, finally we obtain the following sum rules: 

/|M|e Ji^ = dspj/rj{s)e ^ , (13) 



1 
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We carry out operator product expansion to the vacuum condensates adding up 
to dimension-6. In calculation, we take assumption of vacuum saturation for high 
dimension vacuum condensates, they are always factorized to lower condensates with 
vacuum saturation in the QCD sum rules, factorization works well in large N^. limit. 
In this article, we take into account the contributions from the quark condensates, 
mixed condensates, and neglect the contributions from the gluon condensate. In 
calculation, we observe the contributions from the gluon condensate are suppressed 
by large denominators and would not play any significant roles [501 ED ESI [53l El] . 
Furthermore, we neglect the terms proportional to the iriq (= rriu = frid), their 
contributions are of minor importance and can be neglected safely. 



Differentiating the Eq.(13) with respect to -ga, then ehminate the pole residue 
fzi we can obtain two sum rules for the masses of the Z, 



,2 Hl^dspJhis^^ ^ 



Mi = -^ ^ . (15) 

J4m^JspJ/v{s)e ^ 

3 Numerical results and discussions 

The c-quark mass appearing in the perturbative terms (see e.g. Eq.(14)) is usu- 
ally taken to be the pole mass in the QCD sum rules, while the choice of the 
rric in the leading-order coefficients of the higher-dimensional terms is arbitrary 
It is convenient to take the pole mass rric = (1.3 ± 0.1) GeV [S5]- The 



MS mass mc{rnl) relates with the pole mass rh through the relation mc{'m1) 



m 



1 + ^£i^ + {K- 2Cf) (^) + 



, where K depends on the flavor num- 
ber Uf. In this article, we take the approximation rric ~ rh without the ag corrections 
for consistency. The value hsted in the PDG is mciml) = 1.271^^^1 GeV [SZ], it is 
reasonable to take the value rric = (1-3 ± 0.1) GeV in Ref.[55], we also present the 
result with larger uncertainty. The vacuum condensates are scale dependent, the av- 
erage virtuality of the quarks is characterized by the Borel parameter M^, it makes 
sense to choose /^^ = (!?(M^). In this article, the energy scale is taken as /x = 2 GeV, 
(qq) = -(0.26 ± 0.01 GeV)^ {qgsCrGq) = ml{qq), ml = (0.8±0.2) GeV^ 071 ill [56] . 

In the conventional QCD sum rules ^7\ HH], there are two criteria (pole dom- 
inance and convergence of the operator product expansion) for choosing the Borel 
parameter M^ and threshold parameter sq- In calculation, we usually consult the 
experimental data in choosing those parameters. 

The Belle collaboration observed the resonance-like structures Zi and Z2 in the 
TT^Xci invariant mass distribution near 4.1 GeV in the exclusive 5° -^ K'n'^Xci 
decays [IS]. The Breit-Wigner masses and the widths are about Mi = 4051 ± 
14+^0 MeV, Ti = 82+21+47 MeV, M2 = 4248+^^+^f MeV and T2 = 177+^^+^f MeV. 
If they are vector hidden charmed tetraquark states, the central value of the thresh- 
old parameter can be tentatively taken as sq = (4.248 -|- 0.5)^ GeV^ ~ 23GeV2, 
where we choose the separation between the ground states and flrst radial excited 
states to be 0.5 GeV. 

The present experimental knowledge about the phenomenological hadronic spec- 
tral densities of the tetraquark states is rather vague, whether or not there exist 
tetraquark states is not conflrmed with confldence, and no knowledge about the 
high resonances; we can borrow some ideas from the baryon spectra [57j . 

For the octet baryons with the quantum numbers I{J^) = \{.\ ), the mass of the 
proton (the ground state) is Mp = 938 MeV, and the mass of the flrst radial excited 
state A^(1440) (the Roper resonance) is M1440 = (1420-1470) MeV ^ 1440 MeV [57]. 
For the decuplet baryons with the quantum numbers I{J^) = |(| ) , the mass of the 
A(1232) (the ground state) is M1232 = (1231-1233) MeV ^ 1232 MeV, and the mass 



of the first radial excited state A(1600) is Mieoo = (1550 - 1700) MeV ^ 1600 MeV 
[57] . The mass gap between the ground states and first radial excited states can 
be chosen as 0.5 GeV. In this article, the central value of the threshold parameter 
So = 23 GeV^ makes sense. 

However, the threshold parameter sq = 23 GeV^ cannot result in a reasonable 
Borel window, we have to postpone it tentatively to larger values. It is not an 
indication that non-existence of the vector hidden charmed tetraquark states below 
4. 7 GeV; in other words, the QCD sum rules alone cannot indicate (non-) existence 
of the multiquark states strictly. 

If the multiquark states exist indeed, we can release the criterion of pole dom- 
inance and take a more phenomenological analysis with the QCD sum rules. One 
may refuse the value extracted from continuum dominating QCD sum rules as quan- 
titatively reliable if one insists on that the contribution from the pole term should 
be larger than (or about) 50% (for detailed discussions about this subject, one can 
consult Ref.[50]). In the present case, the numerical results indicate that the thresh- 
old parameter sq > 30 GeV^ can lead to possible Borel window, we take the value 
So = (32 ± 1) GeVl 

With the central values of the input parameters, the contribution of the con- 
densate (of the largest dimension) (gg)^ is less than 28% (32%) at the value M^ > 

3.4 GeV^, and the contribution decreases quickly to about 10% (11%) at the value 
M^ = 4.5 GeV^ for the C — 6*7^ (C'Ts — C'T/xTs) type interpolating current. In this 
article, the Borel parameter can be taken as M^ > 3.4 GeV^, we expect the operator 
product expansion is convergent. 

The contribution of the pole term is lager than 50% (49%) at the value M^ < 

4.5 GeV^, and the pole contribution is about (50 — 74)% ((49 — 71)%) at the value 
M^ = (3.4 — 4.5) GeV^ for the C — 6*7^ (C75 — 6*7^75) type interpolating current, 
again we take the central values of the input parameters. If we take into account the 
uncertainty of the threshold parameter, so = (32 ± 1) GeV^, the pole contribution 
is about (48 — 77)% ((45 — 75)%). The Borel parameter can be taken as M^ = 
(3.4 — 4.5) GeV^, where two criteria of the QCD sum rules are full filled [171 HH]. 
From the Figs. 1-2, we can see that the sum rules are not stable enough below 
the value M^ = 3.8 GeV^. In the article, the Borel parameter and the threshold 
parameter are taken as M^ = (3.8 — 4.5) GeV^ and So = (32 ± 1) GeV^, respectively. 

For the tetraquark states consist of light fiavors, if the perturbative terms have 
the main contribution (in the conventional QCD sum rules, the perturbative terms 
are always have the main contribution), we can approximate the spectral density 
with the perturbative term (where the A are some numerical coefficients) |58j . 

BmII^-A s^e'^ds = AM^^ t^e-^dt, (16) 

Jo Jo 

take the pole dominance condition. 



j^t^e-^dt 



>50%, (17) 



and obtain the approximated relation, 

to = -^>4.7. (18) 

The superpositions of different interpolating currents can only change the contribu- 
tions from different terms in the operator product expansion, and improve conver- 
gence, they cannot change the leading behavior of the spectral density p{s) oc s"^ of 
the perturbative term [55] . 

This relation is difficult to satisfy for the light flavor tetraquark states [SOI ED ESI 
[531 El], however, it is not an indication that non-existence of the tetraquark states. 
The hidden charmed and bottomed tetraquark states, and open bottomed tetraquark 
states may satisfy the relation, as they always have larger Borel parameter M^ and 
threshold parameter sq [201 [Ml ESI ED] • Although the relation is derived for the light 
flavor quarks in the massless limit, the c and b are heavy quarks. 

For examples, in Ref.[59], the authors take the X(3872) as hidden charmed 
tetraquark state and calculate its mass with the QCD sum rules, the Borel parameter 
and threshold parameter are taken as M^ = (2.0-2.8) GeV^ and sq = (17-18) GeV^; 
in Ref.|3l], the authors take the Z(4430) as hidden charmed molecular state and 
calculate its mass with M^ = (2.5 - 3.1) GeV^ and sq = (23 - 25) GeV^ In those 
sum rules, the relation in Eq.(18) can be well satisfied. 

In this article, sq/M"^ > 6.8, the relation in Eq.(18) is certainly satisfied. The 
relation can serve as an additional constraint in choosing the Borel parameter and 
threshold parameter, it alone cannot lead to satisfactory results, as a number of 
values of the Borel parameter and threshold parameter satisfy the relation. 

Taking into account all uncertainties of the input parameters, finally we obtain 
the values of the masses and pole residues of the Z, which are shown in Figs. (1-2), 

Mz = (5.12 ± 0.15) GeV, 

fz = (1.31 ± 0.26) X 10-^ GeV; (19) 

Mz = (5.16 ±0.16) GeV, 

fz = (1.25 ± 0.25) X 10"^ GeV, (20) 

for the C — Cjfj, type and C75 — Cjfj^-y^ type interpolating currents respectively. In 
numerical calculations, we observe the uncertainties come from the parameter ttt-q 
are very small, while the uncertainties come from the parameters rric, (qq) and Sq 
are comparable with each other. 

If we take larger uncertainty for the pole mass, rric = (1.3±0.2) GeV, the values of 
the mass change to Mz = (5.12±0.28) GeV and Mz = (5.16±0.30) GeV respectively. 
Furthermore, we vary the energy scale from fi^ = m?^ to iJ? = 4.5 GeV^, the central 
value of the mass Mz changes slowly, less than 0.05 GeV. 

The value ^/so — Mz ~ 0.5 GeV happens to be the energy gap between the 
ground states and first radial excited states of the light baryons [57], the threshold 
parameter sq = (32 ± 1) GeV^ makes sense. 
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Figure 1: The mass Mz with variation of the Borel parameter M^, A for the C—C'jfj^ 
type current and B for C75 — 6*7^75 type current. 



At the energy scale /i = 2 GeV, ^ ^ 0.09 [57], if the perturbative 0{as) 
corrections to the perturbative term are companied with large numerical factors, 
1 + ^(s, nic) — , for example, ^(s, rric) > — ~ 10, the contributions may be large. We 
can make a crude estimation by multiplying the perturbative term with a numeri- 
cal factor, say 1 + ^(s,mc)— = 2, the mass Mz decreases slightly, about 0.1 GeV, 
the pole residue changes remarkably. The main contribution comes from the per- 
turbative term, the large corrections in the numerator and denominator cancel out 
with each other (see Eq.(15)). In fact, the C,{s, rric) are complicated functions of the 
energy s and the mass rric, such a crude estimation may underestimate the 0{as) 
corrections, the uncertainties originate from the 0{as) corrections maybe larger. 

The hidden charmed tetraquark state with the spin-parity 1~ lie in the region 
(5.0 — 5.3) GeV, the mesons Zi, Z2 or Z (about (4.1 — 4.3) GeV) in the vr+Xci 
invariant mass distribution cannot be vector tetraquark states. The spins of the 
Z are not determined yet, they may be scalar hidden charmed states which may 
lie in the region (4.1 — 4.3) GeV j49j, or more likely, they are hadro-charmonium 
resonances [H], P-wave D^D^ molecular states, or D^D^ + D^D\ molecular states; 
more experimental data are still needed to identify them. 



4 Conclusion 

In this article, we assume that there exist hidden charmed tetraquark states with 
with the spin-parity J^ = 1~, and calculate their masses with the QCD sum 
rules. The numerical result indicates that the masses of the vector hidden charmed 
tetraquark states are about Mz = (5.12 ± 0.15) GeV or Mz = (5.16 ± 0.16) GeV, 
which are inconsistent with the experimental data on the tt~^Xc1 invariant mass dis- 
tribution. The hidden charmed mesons Zi, Z^ or Z may be scalar hidden charmed 
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Figure 2: The pole residue fz with variation of the Borel parameter M^, A for the 
C — Cj^ type current and B for C75 — Cj^'j^ type current. 

states, hadro-charmonium resonances or molecular states. 
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